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We  consider  a  DS/CDMA  communication  system  with  K  active  users.  Each  user 
is  assigned  a  unique  spreading  sequence  {ck  (t)}k=1  with  finite  support,  ck  (t)  =  0,  t  </ 
[0,  Tb] .  Tb  is  the  bit  duration.  The  transmitted  signal  is  given  by  s  ( t )  =  ^2k=1  Akbkck  ( t ) , 
t  G  [0 ,7V],  where  Ak  is  the  kth  user’s  amplitude,  bk  is  the  binary  data  bit  transmitted 
by  kth  user.  {bk}k=1  are  assumed  to  be  i.i.d.  Bernoulli  distributed.  Under  BPSK 
modulation,  bk  G  {  —  1,  +1}.  The  size  of  each  vector  is  equal  to  the  length  of  spreading 
sequence  or  processing  gain  L.  Then,  each  vector  can  be  written  as  s  =  Akbkck  = 
CAb,  where  s  is  the  transmitter  signal  vector,  the  dimension  of  s  is  equal  to  L  at  chip 
interval  Tc.  C  G  RLxK ,  C  =  [c1>c2,  •  •  •  , cK\.  CiC2,--  -  , cK  is  the  spreading  sequence 
of  each  user,  respectively.  We  consider  a  MIMO  channel  with  N  transmitted  and  M 
received  antennas.  Assuming  that  the  detector  has  the  knowledge  of  carrier  frequency 
and  phase  shift  and  thus  cancelled  out  in  the  receiver  front  end  and  the  fading  is  slow 
such  that  the  channel  remains  constant  during  one  bit  time  slot  of  J  chip  intervals.  The 
received  signal  at  the  antenna  array  can  be  written  as 

x  (j)  =  Hs  (j)  +n  (j) ,  j  =  1,  •  •  •  ,  J  (1) 

where  H  is  the  M  x  N  MIMO  channel  matrix  coupling  the  transmitter  and  receiver 
antennas,  n  is  white  Gaussian  random  vector.  The  one  bit  received  signal  can  be 
written  as 

X  =  [  x(l)  x  (2)  x(J)  ] 

=  H  [  s  (1)  s  (2)  •••  s(J)]  +  [n(l)  n  (2)  •••  n  ( J)  ]  =  HS  +  N  (2) 

where  N  consists  of  n  (1) ,  •  •  •  ,  n  ( J)  and  is  also  white  Gaussian  random  matrix.  The 
array  output  signal  r  (j )  is  linearly  combined  through  a  complex  weight  matrix  W  to 
yield  the  array  output  signal  that  can  be  written  as  r  (j)  =  WJ/x(j),  j  =  1,  •  •  •  ,  J, 
where  W  is  M  x  N  weight  matrix.  From  (1),  we  reconstitute  the  array  output  signal. 
Buffering  the  array  output  signal  into  blocks  of  length  J,  each  block  can  be  written  in 
vector  form 

r=  W  x  (3) 

where  W  is  JM  x  JN  weight  matrix  and  x  is  JM  x  1  received  signal  vector  that  can 
be  expressed  as 

x=  Hs  +  n  =  HCb  +  n  (4) 

where  H  is  J M  x  .IN  MIMO  channel  matrix. 

TLS  Detector 

The  robust  detector  is  based  on  the  TLS  algorithm[3]  and  we  name  it  as  the  “TLS 
detector [4]”.  As  signal  propagates  through  the  MIMO  channel,  it  is  usually  subject  to 
multipath  fading  such  that  it  arrives  at  the  destination  along  a  number  of  different  paths. 
This  leads  to  the  signature  waveform  distortion [5].  Combine  (3)  and  (4),  the  antenna 
array  output  signal  can  be  written  as 

r  =  W^HCb  +  W^n  (5) 
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The  summary  of  the  TLS-based  algorithm  can  be  described  as  follows,  the  optimum  E,s 
is  composed  of  the  K  largest  (nonzero)  eigenvectors  of  the  matrix  [  C  ?  ]  [  C  r  ]  : 


1.  Perform  eigendecomposition  on 

C  r 

C  f 

to  obtain  Es 

2. Find  projection  operator  Pc  =  E.SE(/ 

3. The  TLS  detector  b TLs  =  (C7  P 

ccy'c 

TPC?/L 

4. The  decision  rule  bTLg  =  sgn  (TPLSr / X) 

ESPRIT  Detector 

The  summary  of  the  ESPRIT-based  algorithm  can  be  described  as  follows: 


1.  Obtain  an  estimate  Rxx  of  R,  ,  from  the  measurements  x 

2.  Perform  eigendecomposit-ion  on  Rxx  to  obtain  Us 
3. Solve  U,,  =  'EUs2  lor  T 


4.LS  solution  4fLS  =  (ufiU^) 
TLS  solution  ^tls  =  — 
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5. Obtain  eigenvalue  Ai,  A2,  •  •  •  ,  Xn  of  'E Lg  and 
b.Find  0r  =  sin^1  (jgj ,  n  =  1,  •  •  •  V 
7.  Obtain  channel  estimates  H  and  H 


8.  The  ESPRIT  detector  b esprit  =  (wHHc)  +  r/h=  Ttespritv/L 

9.  The  decision  rule  b esprit  =  sgn  (T^SPRITr /L) 


LMS  Receiver 


We  propose  an  adaptive  receiver  that  is  based  on  the  LMS  algorithm  and  we  name 
it  as  the  “LMS  receiver” .  The  despreader  output  b  (n)  can  be  written  as  b  (n)  = 
CHr  (n)  /L  =  C2WH  (n)  x  (n)  /L.  A  linear  receiver  that  comprising  a  weight  matrix 

W  operates  on  the  input  vector  x  to  yield  the  output  b.  The  weight  matrix  W  is  chosen 
to  minimize  the  cost  function 


Jlms  =  E 


CHWH 


nix 


(n)  /L  —  b  (n) 


=  -E[l[e(«)H2] 


(6) 


The  adaptive  algorithm  is  as  W  (n  +  1)  =  W  (n)  —  2/ix  (n)  e11  ( n )  C H/L.  In  order  to 
improve  the  LMS-type  receiver,  we  propose  an  advanced  scheme  named  as  the  “LMS 
sequence  revising  (LMS-SR)  receiver” .  The  new  structure  uses  LMS  algorithm  to  update 
despreading  sequence.  Taking  gradient  with  respect  to  the  despreading  sequence,  we 
get  9JqqS  =  2WH  (n)  x  (n)  eH  (n)  /L.  The  sequence  update  equation  is  C  (n  +  1)  = 
C  (n)  - 


ESPRIT  Based  LMS  Receiver 


The  summary  of  the  ESPRIT-LMS  algorithm  can  be  described  as  follows: 
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1 .  Obtain  an  estimate  Rxx  of  R,  ,  from  the  measurements  x 

2. Perform  eigendecomposition  on  Rxx  to  obtain  Us 

3. Solve  Usi  =  Ws2  for  \l> 

4.  Obtain  eigenvalues  Ai,  A2,  •  • 

•  ,  Aat  of  &ls  and  &tls 

5. Find  0,  =  sin"1  ( ln*"  ) 
\  ]kcd  y 

1 ,71  = 

C---N 

6. Obtain  channel  estimates  H  and  H 

7.Initialization  W (0)  = 

(H+) 

\ 

8. The  error  signal  e  =  b 

-b 

9. Weight  update  W  (n  +  1)  = 

W  (n)  -  2/ixe"C H/L 

10. Spreading  sequence  update  C  (n  +  1)  =  C  (n)  —  2/./W H5teH/L 

Simulation  Results 


We  compare  the  performance  between  the  six  proposed  receivers  (TLS-based  detec¬ 
tor,  ESPRIT-based  detector  and  their  adaptive  detectors)  under  the  following  MIMO 
DS/CDMA  environments: 


Transmit  ant.  size 

Receive  ant  size 

User  number 

multipath  number 

case  ‘F 

N  =  4 

M  =  16 

K  =  4 

P  =  5 

case  ‘2’ 

N  =  4 

M  =  16 

K  =  32 

P  =  20 
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Number  of  iteration 


Figure  2:  Bit  error  rate  performance  comparison  between  the  proposed  receivers  in  case 


